Terahertz spectroscopy of spin waves in multiferroic BiFeOs in high magnetic fields 
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We have studied the magnetic field dependence of far-infrared active magnetic modes in a single 
ferroelectric domain BiFeOs crystal at low temperature. The modes soften close to the critical field of 
18.8 T along the [001] (pseudocubic) axis, where the cycloidal structure changes to the homogeneous 
canted antiferromagnetic state and a new strong mode with linear field dependence appears that 
persists at least up to 31 T. This allows us to assign the low field modes to the excitations of the 
cycloid and the high field mode to the antiferromagnetic resonance. We found that the zero field 
spectra at low temperature are different before and after application of high magnetic field and this 
change persists at low temperature. We interpret this as a change from a state with three magnetic 
domains into a state with one domain with the lowest energy in high magnetic field. 

PACS numbers: 75.85. +t, 76.50.+g, 78.30.-j 
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FIG. 1. (color online) Pseudocubic unit cell of BiFeOs show- 
ing the positions of Bi and Fe ions, the ferroelectric polariza- 
tion P, three equivalent directions of the cycloidal ordering 
vector qi, the applied static magnetic field Bo || [001], and 
the wave vector of incident light k together with the electric 
field (e) component of light in two orthogonal polarizations 
that were used in the experiment. 



Multiferroics are materials with more than one fer- 
roic order parameter and they have been in the focus 
of recent research [THZ]. Among multiferroics BiFeOa 
is the only known one that exhibits both ferroelectric 
and magnetic orders already at room temperature. Be- 
low Tcuric ~ HOOK it has a rhombohedrally distorted 
perovskite structure [H], space group R3c, were the iron 
(Fe'^+) and bismuth (Bi^"*") cations are displaced from 
their centrosymmetric positions and produce a polar dis- 
tortion, resulting in a spontaneous ferroelectric polariza- 
tion P along the [111] direction of the pseudocubic unit 
cell [9]. 

BiFeOs shows G-type antiferromagnetic (AFM) order 
below Tjv ~ 640 K in which the Fe'^"'" magnetic moments 
are ordered antiferromagnetically between adjacent pseu- 
docubic (111) planes [lOl [11]. Ferromagnetic ordering 
within the planes is broken by a long-wavelength cy- 
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FIG. 2. (color online) Absorbance spectra of cycloid modes 
in zero field and in e || [110] polarization (a) and in e [| 
[110] polarization (b). Solid curves, ap, were measured after 
applying the high field Bq > 21 T. Dotted curves are initial 
absorbance spectra ai of the zero field cooled sample. 



cloidal spin order [TT] with three possible directions of 
the ordering vector q, || {[1, 1, 0], [0, 1, 1], [1, 0, 1]} _L P, 
see Fig.[l] The spins of a cycloid are in the plane 
determined by P and q^. 

The cycloidal order in BiFeOa is induced by a weak 
relativistic Dzyaloshinskii-Moriya interaction with cou- 
pling constant D^, _L P and Dj, _L q, [T^. Another 
Dzyaloshinskii-Moriya interaction D || [111] cants spins 
out of the cycloid plane [5J [H] and the ferromagnetic or- 
dering of this canted moment has been verified by a neu- 
tron scattering experiment [13. . High resolution neutron 
scattering shows that the magnetic ground state ordering 
in BiFeOa does not change in zero field on cooling from 
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TABLE I. Fitted positions and areas of absorption lines from Fig. [2] Ai and Af are the line areas before and after application 
of Bo > 21 T. Fit errors are three times standard error. Mode labels are according to Fishman et al. [27]. The $_i mode 
position is an extrapolated value of the magnetic field dependence to OT, Figlll Aq^ and Aq^ ^ 3 s-re theoretical line areas in 
arbitrary units assuming a magnetic monodomain with qi and an equal popuTation of qi, q2 and qa domains, respectively. 
The line areas of three modes are *_i : ^ = J2^rl^\{q^\M,^\0)\'^hl^, ^-+1 : ^ = Ei '^41 (<?i|^-, |0)p/i^, + \{q^\PyM\'^4,]' 
$±2 : A = Wi|{(7i|MaJO)p/iy., where Xi || q^, Zi \\ [111] and = Zi x Xi, populations m + n2 + m = I and |hj = |e| = 1. 
Matrix elements are |27] \{qi\M^,\0)\ = 4,\{qi\My^\0)\ = 0.8, |(<ji|M^JO)| = 3.6 at K = 0.002meV, p = \{q^\Py,\0)\. 
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300K to 4K [HHS]. 

Single ion anisotropy along the easy axis [111] intro- 
duces anharmonicity ^121 117] , but in zero magnetic field 
the cycloid is only slightly anharmonic [121 HI] • Exter- 
nal magnetic field contributes to the effective single ion 
anisotropy (19. and induces a metamagnetic transition 
[12] at a critical field B^. « 19 T, where the cycloidal 
order changes to a coUinear AFM spin order [20] • The 
unwinding of the cycloid reduces the electric polarization 
[21j and creates a small macroscopic spontaneous mag- 
netization induced bv D || [111] [1^1 [2T]. 

In BiFeOa it has been shown that electric field flips 
polarization P together with the cycloidal plane and 
changes the populations of magnetic domains P^H^ . To 
our best knowledge there is no information whether the 
applied magnetic field reorients magnetic domains in a 
single crystals or not, although the simultaneous action 
of electric and magnetic, internal or external fields, is 
crucial in the operation of exchange-biased devices |2S] . 

The eigenfrequencies of magnetic modes are sensi- 
tive to anisotropic magnetic interactions and these in- 
teractions are important to understand the microscopic 
models behind magneto-electric coupling in multiferroics. 
The eigenspectrum of BiFeOs cycloids was calculated by 
de Sousa and Moore and with the addition of single- 
ion easy-axis anisotropy by Fishman et al. |27j in T and 
in applied electric field by Rovillain et al.[35]. Spec- 
troscopic techniques what measure the eigenspectrum 
of magnetic modes, especially if they can be combined 
with external fields what compete with internal fields, are 
valuable tools. Most of the INS US] [3D], Raman [5T|[5^ 
and THz [33J [M] spectroscopy studies on BiFeOs were in 
zero applied field. The Raman work demonstrated that 
the Raman-active magnon frequencies depend strongly 
on applied electric field [28|- The high field ESR was 
done in magnetic fields up to 25 T, but was limited to 
frequencies lower than the main cycloid modes and one 
of the AFM modes. 

In this work, we present THz spectra of BiFeOa single 



crystals at low temperature and Bq < 31 T which allow 
us to identify the excitations of the cycloid and follow 
their magnetic field dependence until the cycloidal order 
is destroyed in high magnetic field. We show that mag- 
netic field changes the distribution of magnetic domains. 

The (001) face single crystal BiFeOa sample was grown 
using a Bi203 flux |35| . It has a thickness of 0.37mm and 
it contains a single ferroelectric domain along the [111] 
axis, Fig.[l] checked by an optical rectification experi- 
ment |36| . 

The sample was zero field cooled and spectra were mea- 
sured in Faraday configuration with the magnetic field 
along the [001] axis. Up to 12 T spectra were measured 
at 4K in Tallinn with a Martin-Puplett spectrometer 
and a 0.3 K bolometer [37' using a spectral resolution 
of 0.2 cm~^. Spectra from 12 T up to 31 T were mea- 
sured in Nijmegen High Field Magnet Laboratory at 2 K 
using a Bruker IFS 113v spectrometer and a 1.6 K sili- 
con bolometer and spectral resolution of 0.43 cm^^; the 
spectra were averaged for 15 minutes at each field. There 
was a linear polarizer in front of the sample to control the 
polarization of light. 

We measured absorbance spectra in magnetic field 
with the reference spectrum in zero field. This method 
gave excellent spectra of magnetic field dependent lines. 
From the relative absorbance spectra in fields above 21 T 
(after applying 30 T) we extracted the zero field absorp- 
tion lines, solid curves in Fig.[2]and fitted them. The zero 
field fit results, parameters in Tabled were added to the 
measured relative spectra in magnetic fields. The result, 
absolute absorbance spectra in fields, are shown in Fig.[3j 
and fitted line positions and areas in Fig.|4] 

We observed a change in the zero field spectra after 
applying high field at low temperature, Fig.[2j The zero 
field spectrum stayed the same after applying high field 
in the opposite direction. The initial zero field line inten- 
sities «!, measured on the zero field cooled sample, were 
recovered after warming the sample to 300 K. This is an 
evidence that different magnetic domains exist. 
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FIG. 3. (color online) Magnetic field dependence of ab- 
sorbance spectra of BiFeOa in e || [110] (black) and e || [110] 
(red) polarizations. Zero field spectra are fits as described in 
the text. The inset shows a detailed view in the spectral range 
of the AFM resonance close to the metamagnetic transition, 
Be = 18.8 T. 



FIG. 4. (color online) IVIagnetic field dependence of magnon 
modes in the THz spectrum of BiFeOa at low temperature. 
The areas of triangles and circles are proportional to the ab- 
sorption line areas. Vertical dashed lines mark the metam- 
agnetic transition at Be = 18.8 T and another critical field 
Ba ~ 5.5 T. The solid lines are the fit of our data and ESR 
data [38| (squares) above 19 T to a model from Ref.l38l 



When the cycloidal order is destroyed above the critical 
field Be a collincar AFM order is established where the 
spins are in the (111) plane [20] with moments Mi and 
M2 antiparallel in neighboring planes and with a small 
canted moment in the (111) plane, M = Mi -I- M2, in- 
duced by D II [111] 13 m]. The energy is at the minimum 
when M || Bq and for the field [0,0, Bq] this happens for 
the qi domain and the AFIVI vector L = Mi — M2 will 
be aligned along qi. Even if the high field is removed, 
the selected q will remain and a monodomain cycloidal 
state is created. 

We found that the change in the zero field spectra oc- 
curs already by 12 T and higher fields that destroy the 
cycloid do not change the zero field lines any more. Ap- 
parently the qi cycloid where the field is mostly per- 
pendicular to the cycloid plane has a lower energy than 
the cycloids q2 and qs where the magnetic field is closer 



to the cycloid planes. This means that a net magnetic 
moment perpendicular to the cycloid plane is created by 
applied field. Indeed, a non-linear dependence of M{Bq) 
has been observed below 5 T if the direction of Bo devi- 
ates from the [111] direction [21]. The results shown in 
Fig.|3] and [4] were obtained from spectra measured after 
the sample had been in high field > 12T and in further 
discussion we assume that qi is the only cycloid present. 

The zero field line parameters are compared to the cal- 
culation of Fishman et al. [17] in Table|l] using their no- 
tation for the cycloid modes. The largest absorption line 
area is that of the $±2 hue in e ]] [110] polarization and 
it is quite strong also in the other polarization, contrary 
to the theory what predicts the smallest matrix element 
among magnetic dipole transitions for this mode. The 
polarization dependence of the ^'_i and 5'+i modes is 
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better described by the calculation. 

The spectra in field and a detailed field dependence of 
mode frequencies and areas are presented in Fig.|3]and|4] 
We see that more modes, too weak to be detected in zero 
field, coincide with modes predicted in Ref.[?fl llcm~^ 
with $_i and 32cm~^ with ^±2 mode at g = 2, and the 
two modes between 40 and 45 cm"-'^ are even higher order 
excitations. 

The three main modes, ^'^i, ^+1 and $±2 change only 
slightly with increasing magnetic field until a critical field 
Ba ~ 5.5 T is reached, where a discontinuity on several 
modes and a smooth change in the slope of the mode 
is observed, Fig. [4] These changes might be associated 
with the change in the magnetic structure. Indeed, the 
magnetization measurements on similar samples show a 
non-linear behavior of M{Bq) below 5T and almost lin- 
ear behavior above [H]. However, the THz spectra do 
not show any anomalies at 10 T seen by optical measure- 
ments [39] . 

Modes soften before reaching the metamagnetic tran- 
sition at Be, except ^P+i which seems to merge with the 
softening $±2 at about 18 T. There is an intriguing pos- 
sibility that close to transition between 18.6 and 18.8 T, 
see inset to Fig. [3] and |4| one of the cycloid resonances, 
^+1 or $±2 as labeled in zero field, coexists with the 
AFM resonance. This means that in a narrow field in- 
terval the spin structure supports both, cycloid-like and 
AFM-like mode. The coexistence of two phases is ruled 
out since the metamagnetic transition in BiFeOs is not 
the first order phase transition nor similar to a spin flop 
transition in ordinary antiferromagnets }19j and we did 
not observe any hysteresis effects between 18 and 19 T 
as reported earlier |38| . In THz spectra there is only 
one resonance line above 18.8 T and we assign this value 
to the critical field Be of the metamagnetic transition in 
BiFeOg at 2K and Bq || [001]. 

From the ESR [38] and THz spectra we have both 
AFM resonances in hand and we fit them to the model 
described in [35]. The magnetic field Bq projected on 
the axes of the qi cycloid is Bo[0, \/2/3, \/l/3]. We 
get the following parameters: gyromagnetic ratio 7 = 
(1.72±0.01) X 10''rad(sOe)-\ Ku/x± = (1.06±0.02) x 
XQio erg/cm^, H^m = (105 ± 2)kOe, see the solid lines 
in Fig.|4] Ku is the uniaxial magnetic anisotropy con- 
stant and X-L is the susceptibility perpendicular to the 
L. The canted moment x^H^m = Mi -I- M2 is in- 
duced by D. We use x± = 3.1 x 10~^ from the high 
field magnetization measurement [21] and get for the 
magnetic anisotropy constant (energy density) = 
3.3 X 10^ erg cm""^. We get for the same quantity 
{KSD ~ 3.2 X 10^ erg cm"'^ from the microscopic model 
[27] where K = 0.002 me V and assuming 5*^ = 5/2. The 
canted moment estimated from the AFM resonance spec- 
tra is XJ-Hdm — 0.022/iB/Fe. This should be compared 
to 0.05/xb derived from the interpolation of the high 
field magnetization to the zero field (0.04/iB/Fe projected 



on [0, ^1/3] direction) [5T] and to averaged local 

magnetization 0.06^B/Fe in the cycloidal phase 

In summary, the cycloid modes soften close to the 
metamagnetic transition at Be — 18.8 T and are replaced 
by two AFM modes, one observed earlier by ESR [38j . 
The parameters derived from the AFM resonances are in 
reasonable agreement with neutron scattering [16] and 
magnetization |21| measurements and with the micro- 
scopic spin model [27) . This model describes properly 
the resonance frequencies of the cycloid modes in zero 
magnetic field. However, the assumption that these tran- 
sitions are only magnetic dipole active is not sufficient to 
account for the observed polarization dependence. The- 
ory is needed to describe the Bq dependence of the cy- 
cloid excitations and the observed discontinuity of mode 
dispersion in magnetic field at 5.5 T calls for the need to 
revise the minimal magnetic model in BiFeOs. Our work 
demonstrates that in addition to electric field [23, the 
control of magnetic domains with magnetic field is pos- 
sible; a neutron scattering experiment where both fields 
are applied is needed to understand multiferroic proper- 
ties of BiFeOa. 
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